Parametric Ward-Takahashi identity in disordered systems and the integral identity 
associated with the Calogero-Sutherland model 
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By utilizing the symmetric property known as the Ward-Takahashi identity in disordered sys- 
tems, we explore the novel symmetry relations which hold in one-dimensional systems with inverse 
square interaction (the Calogero-Sutherland model). The identities emerge totally from the algebraic 
structure of the model. They show that the dynamical correlators are connected with one another, 
involving the higher-order integrals of motion. We obtain the result for the coupling strengths 
A = 1/2, 1, and 2, and conjecture that a similar relation may hold for arbitrary rational A. 
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I. INTRODUCTION 

The Calogero-Sutherland model (CSM) describes N 
fermions located on a ring of the perimeter L with pair- 
wise inverse square interactions 0B]. The Hamiltonian 
is given by 



H, 



CSM 



1 v^ , n^A(A-l) 
— / Pi ^ ^ 



2m 



E 



sm 



(1) 



where pi = —ihd/dri, rij — ri — tj, and (f> — tt/L. 
(The usual convention Ti = I and m — 1/2 is adopted 
hereafter.) The exact ground state wavefunction of the 
model is given by a Jastrow form ni<7^in- (^r^), and 
excited states are known to be expressed in terms of the 
symmetric polynomials called the Jack polynomials [gj. 
While there is a long history of exactly solvable mod- 
els in one-dimensional many-body systems, this model is 
the only family known so far where its dynamical correla- 
tion functions can be evaluated exactly. The dynamical 
density-density correlator (p(r, t)p(0, 0)) in the thermo- 
dynamic limit was evaluated analytically for integer and 
rational values of A = p/q ^Q- A striking simplicity of 
the result emerges after taking the thermodynamic limit, 
which was attained through a lot of mathematical effort. 
The Calogero-Sutherland model is closely connected 
with the random matrix theory (RMT), which has suc- 
cessfully been applied to describe the universal charac- 
teristics in quantum chaotic systems, such as compound 
nuclei, quantum billiards and quantum dots. In the ther- 
modynamic limit, the Jastrow form of the wave function 
immediately enables the ground state average to be iden- 
tified with the average over Wigner-Dyson ensembles of 
random matrices for coupling strengths A = f3/2 — 1/2 
(orthogonal), 1 (unitary), and 2 (symplectic). The spec- 
tral correlator was generalized to account for spectra that 
disperse as a function of some external tunable param- 
eter. Surprisingly, it was found that this parametric 
two-level correlator is identical to the dynamical density- 
density correlator of CSM ^^ ■ Though this "mapping" 
between RMT and CSM is available only for the three 
special values of the coupling strengths, it can serve as a 



rich source of various useful insights in CSM. (See, e.g., 
Rcf. Q for recent work in this direction.) 

In this paper, we utilize the mapping to examine the 
symmetric relation associated with the dynamical corre- 
lations in CSM. To do so, we extend the Ward-Takahashi 
identity in disordered systems to incorporate the para- 
metric correlations with the help of the supermatrix 
method. In contrast to the Jack polynomial technique, 
this method is suitable to investigate physical quantities 
in the bulk limit. It also transparently provides symmet- 
ric properties to dynamical correlation functions, which 
we will explore. 



II. WARD-TAKAHASHI IDENTITY 

Our starting point is the Ward-Takahashi identity in 
quantum dots or RMT. It asserts that, for the retarded 



and advanced Green functions G 
the identity 



R,A 



(E - H ±ioy 



Tr 






2ni 



A{Ei ~ E2 



(2) 



be satisfied in disordered systems, where A is the mean 
level spacing and ~^^^ denotes the averaging over impu- 
rity configurations or random matrices. Eq. (0) results 
from the unitarity of the system, so it should be possi- 
ble to extend this identity to incorporate the paramet- 
ric dependence. Although Eq. (||) itself can be proved 
straightforwardly by inserting the complete diagonalized 
basis between G^ and G^, such a route of derivation 
is no longer achieved when they carry different external 
parameters, since we can make no common diagonalized 
basis. To extend Eq. (||) to such situation, we should 
take account of the unitarity of the system explicitly. To 
do so, we resort to the supermatrix method |10(] which 
translates the unitarity of the system into the symme- 
try of the advanced and retarded components. When we 
present the results in the context of CSM (Eqs. (25 28) 
below), it will be found that the derived identities take 
quite simple forms, but still give nontrivial relations even 
for the free fermion case (A = 1). 



III. SUMMARY OF THE CORRELATORS IN 
RMT AND CSM 

A. Parametric correlators of RMT 

To make our discussion concrete, take the Hamiltonian 

H{X) = Ho+ X$, (3) 

where Hq is a random matrix belonging to one of the 
Dyson ensembles, and <I> is a fixed traceless member of 
the same ensemble. By use of the retarded and advanced 
Green functions 



G'yxir,r') = {r\{E-HiX)±iO)-'y) 



(4) 



we define the following two kinds of universal correlation 
functions k{uj,x) and n{uj,x) ||Tl|: 

Hlj,x) = -^ + ^ Jdrdr' G^ir,r)Gi{r',r') (5) 

n{uj, x) = —^ drdr' Gf (r, r')G^{r', r) (6) 

2n'' J 



where the suffix i — 1,2 denotes {Ei, Xi) and s is the level 
degeneracy which takes account of the Kramers doublets 
for the symplectic case. Rescaled parameters for the en- 
ergy and the external parameter were introduced to re- 
veal the universality by 



w = (Si - E2)/A; 
C(0) = Z ~ 



x"" = GiO)iX, - X2)^ 



idE„{X)/dXY 



For the orthogonal, unitary, and symplectic ensembles, 
the analytical answers for k(uj,x) and n{uj,x) have been 
obtained [|ll|Jl2| ]. To present the results simultaneously 
for all three ensembles, the integral variables introduced 
in Ref. ||] are convenient. By assigning A — p/q = 1/2 
for the orthogonal, A = p/q = 1 for the unitary, and 
-^ — p/q = 2 for the symplectic symmetries (p and q are 
coprimes), they are presented by 



k{uj,x) = X 


' q2 ^iQ^-Ex'^ 11 


J 


(7) 


n{LU, x) = X 


'^^iQw-Ex'^/2 




(8) 


The integration X[- ■ ■] is defined by 




1 i-oo 

X[--]^G\\ / dxA 
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Xr^Vj}) (• 


•), (9) 



J=l' 



g = (27r) 



E = (271) 



i=i j=i 






(10a) 
(10b) 



and epix) = x{x + A) and tH{y) = Aj/(1 — y). The 
numerical constant G and the form factor F(X\{xi,yj}) 
were given by 



G 



x2p{q-i) r2(p)r'?(A)rp(i) 



2^ Vein r'b - A(* - i))n r2(i - ^) 

i=l 3=1 



-, (11) 



FiX\{x,,yj}y 



xf[ep{x.r-'f[eHiy,y/'-'. (12) 

i=l j=l 



In the supermatrix formulation, F(X\{xi,yj}) emerges 
as a Jacobian for the integration which is completely de- 
termined from the structure of the underlying graded- 
symmetric space. 



B. Connection with CSM 

The direct connection between the parametric corre- 
lations of RMT and dynamical correlations of CSM is 
provided when we substitute uj —> r and x'^ /2 —> t (r is 
the Euclidean time) p-Q|. When we make this replace- 
ment in A;(w, x), it immediately reproduces the dynamical 
density-density correlator (p(r, r)/9(0,0)) for A = 1/2, 1, 
and 2, i.e., 



(p(r, t)p(0, 0)) = J [g2 cos(gr) e"^" 



(13) 



The other function n(w,x) is found to be related to the 
dynamical current-current correlator of CSM M , 



{Jir,T)JiO,0))=X[E^cos{Qr)e- 



-EtI 



(14) 



Since the Ward-Takahashi identity Eq. (0) states 
X[E e'^'^'^] — —l/{iTruj), it characterizes the current- 
current correlator of CSM rather than the density-density 
correlator. 



IV. DERIVATIONS AND RESULTS 

Now we present how we can extend and derive the 
Ward-Takahashi identity to the case for finite x, or dy- 
namical correlations. To avoid the notational confusion, 
we use (w, x) of RMT, instead (r, r) of CSM, but by 
substituting ui ^ r and x^/2 -^ r, we can obtain the 
corresponding expressions for CSM on each step. We fol- 
low Refs. [|^,0 to derive the Ward-Takahashi identity 
within the framework of the supermatrix method. The 
basic underlying idea is to translate the unitarity of the 



system into the hyperbolic symmetry between the ad- 
vanced and retarded components. (See Eq. (^0|) below.) 
Consider the generating function for k(uj,x) and 
n{uj,x)hi the supermatrix nonlincar-cr model formula- 
tion 0, 



Zj = (exp[STr(QJ)])g, 



(15) 



where Q is the 8x8 supermatrix satisfying Q^ = 1 and 
its explicit structure of Q can be found in Ref. ||l^ . The 
supertrace STr is defined by STr(- • •) = Tr[{kF-kB){- ■ •)] 
where ka is a projector either onto the Bose space {a = 
B) or onto the Fermi space {a = F). The source matrix 
J is chosen as (a and h are c-numbers) 



J ^ {a + bT,i)Aka. 



(16) 



Note that we are allowed to use either one to generate 
k{ijj, x) and n(w, x). A and Si are 8x8 matrices defined 
by 



A 



U 

-u 



Si 



I4 
I4 



(17) 



The average (• ■ •)g denotes the integral /_DQ (• • •) e^^I^l. 
Corresponding to A = p/q = 1/2 (orthogonal), 1 (uni- 
tary), and 2 (symplectic), F[Q] is equal to 



F[Q]=p ^!^STr(QA)- 



16 



-STr(QA)^ 



(18) 



When we make the infinitesimal rotation on the saddle- 
point manifold as Q ^ Q' = (1 - 5T)Q{1 + ST), Zj re- 
mains invariant because of the integration over the man- 
ifold Q. Hence 



5Zj = ( eST'S-' STr 



2 2 9 



[J,Q]-^[A,Q] 



{(AQ)2-(QA)2} Ut 



= 0. 



(19) 



Although this identity holds for arbitrary infinitesimal 
rotations 6T, we particularly choose (a = _B or F) 



ST oc fcctSi. 



(20) 



The choice reflects [/(1, 1) symmetry of the advanced and 
retarded components within the bosonic or fermionic sec- 
tor. This hyperbolic symmetry is responsible for pro- 
ducing the Ward-Takahashi identity in the supermatrix 
method. After substituting Eq. ( po[ ) for ST, we have 



bq, + l^i^i^,, + ^,3 ) e-^+^- 



0, 



(21) 



where we define 



qi = STT[k^AQ], 

92 = STr [A:,EiAQ] , 
<73 = STr[fc„Si(AQ)2] 



(22a) 
(22b) 
(22c) 



Note that the correlator k(LL},x) and n{uj,x) are related 

by 



1 

16 
1 

16 



{{12?)q 



1 + k{ijj, x) 
±n{uj, x) 



(23a) 
(23b) 



Depending on the choice of ks or kp, we have positive 
or negative sign in front of n{uj,x). 

From Eq. (|2l|), we can readily derive a sequence of 
integral identities by comparing each coefficient of poly- 
nomials of a and b. Not all of them, however, produce 
nontrivial integral identities. We can show that the co- 
efficients of a'^b'^ and a^b'^ vanish trivially. The first non- 
trivial identity comes from the coefficient of a°6^, i.e.. 



1 
P 



{qi)c 






(12? 



X^n^x^ 



-92 93 



(24) 



After some straightforward but rather lengthy evaluation 
of the supermatrix integration for all three values of A, 
we obtain the result which can be summarized as follows 
(restoring w — > r and x'^ /2 — > t): 



1 



n 

where we define 
9 



-=j[{-irE + Tl3): 



r-ET 



(25) 



(27r)" 



^a;,(xi-|-A)(2a;, -l-A) 



n-2 



.j=l 






(26) 



Eq. (25) serves as the extension of the Ward-Takahashi 
identity Eq. (0), and consists of the main result of the 
paper. 

We can go on to the higher-order identity from 
Eq. (pi|), on principle, but the evaluation of the integra- 
tion becomes harder and harder to complete. Among the 
second-order polynomials of a and b, we can confirm that 
only the coefficient of a^b^ gives the nontrivial relation: 



-(ill)' +(12)'), 



— 91(92) 



X^TT^X^ 



-919293 



Q 

(27) 



However, as we see from Eq. (23b), this will depend on 
the values of A {p and q) as well as the choice of a = i3 or 
F. Hence for each value of A, we have two integral iden- 
tities. From these, we can seek an interesting form of the 



identity which seems the direct extension of Eq. (25) 
which can be presented by 



1+X 



(Q 



2 , 1 P^\ ^iQr^Er 



— l[{~-trh + TL 



P 



4je 



„iUr-l^T 



(28) 



Note that the coefficients of Eq. (E^) are deduced to re- 
produce the actual results of A = 1/2, 1, and 2. 



V. DISCUSSION 



There are known multiple integral identities which are 
associated with CSM. They are called the Selberg inte- 
grals |l5|| , and their generalization by Dotsenko and Fa- 
teev [ |16| are particularly useful. For instance, they were 
used to determine the correct normalization factor of 
the correlation functions ||l^. We also mention that the 
Dotsenko-Fateev integral can provide a systematic means 
to evaluate a certain correlation function which showed 
up in disordered systems p8| . However, we emphasize 
that those integral formulae are not powerful enough 
to explain Eq. (Eq), because they can be applied only 
when the integrands are polynomials. The simple form 
of the identity Eq. (|2^) may suggest that these known 
multiple integral formulae be extended somehow for the 
case involving an exponential factor such as e^Q^-^"^, 
We remark that the derived integral identity Eq. ( |25| ) 
is not trivial at all from the mathematical point of view, 
even for the simplest case of the free fermion (A = 1), 
though we can convince ourselves of its correctness, e.g., 
by checking the asymptotics, or evaluating for small r 
expansion. 

The quantities /„ (for n > 3) correspond to the higher- 
order integrals of motion of CSM, as well as Q and 
I2 = E. To see this transparently, identify the veloci- 
ties Vi {vj ) for particles (holes) by 



Vj ^Vsil-2yj), 



(29a) 
(29b) 



where Vg = TrXTipo/m = 2ttX is the sound velocity [[l9[ . 
Since the velocity (rapidity) is the conserved quantity of 
CSM, 



9 

i=l 



P 



(30) 



should act as the integrals of motion, so does /„ = 
J„/2-27r2AV„_2- In Ref. @, a few kinds of the higher- 
order integrals of motion were investigated. Although the 
similarity of their look, the direct connection with /„ in 
Eq. (Eq) is missing at present. 



VI. CONCLUSION 

In conclusion, we have derived the Ward-Takahashi 
identity for the parametric correlations of RMT. By do- 
ing so, it was shown that there exist novel integral iden- 
tities which are associated with the dynamical correla- 
tions of CSM. It is remarked that they amount to a new 
generalization of the Selberg integration. As was seen 
from the derivation in the context of RMT, this is the 
manifestation of the unitarity of the system, i.e., the hy- 
perbolic symmetry of the advanced and retarded compo- 
nents. However, its nature and implication in CSM is 
not so clear at present. Since our arguments rely heavily 
on the mapping between RMT and CSM, we can make 
no decisive statement on the validity of the derived inte- 
gral identities for arbitrary rational values of A. We can, 
however, suggest two possible scenarios: Eqs. (p3,E8) are 
(1) true only for A ~ 1/2, 1, 2, or (2) true for all rational 
values of A. If the latter were true, it would remain as a 



future challenge how the integral identities Eqs. (25 



can be deduced from the Jack polynomials, or the W- 
algebra which is known as the symmetry of CSM |2^| . 
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